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Le precede consiste a estimer (9) la mi 

signaux'recus du transmetteur sur un ensemble de n cap- 
teurs de reception, 

a calculer (10) la matrice de correlation inverse R^, a calculer (11) des 
vecteurs d'intercorrelation f^j entre les signaux repus sur l ensemble des n 

capteurs et un signal de replique d(k) connu, a calculer un 
critere de synchronisation (15) en effectuant les produits 
des vecteurs d'intercorrelation 

transposes conjjgues et dc 



le des vecteurs d'intercorrelation 



sr (16) la 



determinee pour placer la synchronisation sur I'echantillon 
pour lequel la valeur des criteres depasse la valeur de seuil 
T|. Pour estimer le decalage Doppler, la sequence de syn- 
chronisation est decoupee en tranches d'un nombre deter- 
mine de symboles pour effectuer des calculs de transfor- 
mation de Fourier rapide sur des vecteurs d'intercorrelation 
determines sur chaque tranche. 
Application: Transmission radio de donnees. 




IIIIIIIIIIIIHWUI 



ion, PLLC - http://www.sughrije.cerr: 



2715488 



1 



Procede et dispositif permettant a un modem de se synchroniser 
sur un transmetteur de donnSes numeriques par voie 
hertzienne en presence de brouilleurs 

5 

La presente invention conceme un procede et un dispositif permet- 
tant a un modem de se synchroniser sur un transmetteur de donnees nume- 
riques par voie hertzienne notamment en presence de brouilleurs. 

Dans les systemes de transmissions de donnees numeriques par voie 

10 hertzienne la synchronisation est effectuee, a partir de sequences de sym- 
boles connus du recepteur et intercalees entre les symboles d'information, 
par correlation entre le signal recu et un signal de reference construit a partir 
des symboles connus. Le resultat est compare a un seuil pour decider de la 
synchronisation. Des systemes permettant la synchronisation d'un modem a 

15 partir d'un seul capteur et d'un seul recepteur sont par exemple connus de 
['article D.V. Sarwate et M.B. Pursely intitule "Cross correlation properties of 
pseudo random and related sequences" paru dans IEEE Vol. 66, pp. 593- 
619, May 1980. Mais les techniques mises en oeuvre qui permettent en 
outre une estimation des caracteristiques du canal de transmission devien- 

20 nent rapidement inefficaces en presence de brouillage. 

Lorsque le rapport signal a bruit n'est pas suffisant pour la prise de 
synchronisation par une technique mono-capteur, des techniques visant a 
accroitre le rapport signal a bruit et a permettre la synchronisation sont utili- 
sees. 

25 Ces techniques peuvent etre regroupees en deux families. Une pre- 

miere famille vise a pretraiter le signal recu avant d'effectuer la synchronisa- 
tion du signal pretraite. La seconde famille vise a ameliorer le rapport signal 
a bruit et a effectuer la synchronisation conjointement. 

Dans les techniques relevant de la premiere famille il faut egalement 

30 distinguer deux categories : 

- Une premiere categorie met en oeuvre un pretraitement mono-cap- 
teur de type excision de brouilleurs avant d'effectuer la synchronisation. Ce 
traitement est connu par exemple de I'article de L.B. MILSTEIN "Interference 
rejection Techniques in spread spectrum communications", paru dans Proc. 

35 IEEE, Vol. 76, pp. 657-661, June 1988. Suivant cet article on cherche a 
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predire la contribution du brouilleur, dont la bande est supposee plus petite 
que la bande du signal, sans predire celle du signal utile. Cela permet 
d'ameliorer les performances de la synchronisation en presence d' interfe- 
rences bande etroite (type FSK). Mais le traitement reste inefficace en pre- 
5 sence d'un brouillage occupant toute la bande utilisee par le signal utile. 
Cela est en partie due au fait que les techniques de pretraitement mono- 
capteur ne permettent un traitement des signaux que dans un espace limite 
a I'axe des temps et a I'axe des frequences. 

- La seconde categorie concerne les techniques mettant en oeuvre un 

10 pretraitement multi-capteurs de rejection, precedant la prise de synchronisa- 
tion. L'idee consiste a mettre en oeuvre un filtrage spatial, discriminant le 
signal et les brouilleurs par leur direction d'arrivee, le temps, le spectre ou la 
puissance comme cela est decrit dans I'ouvrage de R.A. Monzingo, T.W. 
Miller, intitule "Introduction to adaptive arrays", edite par John Wiley and 

15 Sons, New-York, 1980 et la these de P. Chevalier, intitulee "Antenne adap- 
tative : d'une structure lineaire a une structure non lineaire de Volterra", 
These de doctorat, Universite de Paris Sud, Juin 1991 puis de correler le 
signal de reference avec la sortie de Pantenne. Les techniques de filtrage 
spatial discriminant le signal et les brouilleurs par le temps, le spectre ou la 

20 puissance sont generalement sous-optimales, alors que celles exploitant la 
connaissance de la direction d'arrivee du signal necessitent la mise en 
oeuvre d'une technique de goniometrie, couteuse et delicate dans la gamme 
des frequences HF notamment du fait des multi-trajets. 

Pour ces raisons, les techniques de la seconde famille, visant a effec- 

25 tuer la rejection de brouilleurs conjointement avec la prise de synchronisa- 
tion, a partir de la connaissance du signal de reference, apparaissent plus 
prometteuses. Une telle technique a par exemple ete proposee par 
M. Compton dans son article "An adaptive array in a spread spectrum com- 
munication system" paru dans Proc. IEEE, vol. 66, pp. 289-298, March 1978, 

30 pour la synchronisation de signaux PN, a partir d'une architecture analogi- 
que integrant des "Delay Lock Loop". 

Une seconde technique fut proposee par L.E. Brennan et I.S. Reed 
dans leur article intitule "An adaptive array signal processing algorithm for 
communications", paru dans IEEE Trans-Aero-Syst. vol. AES-18, n°1, pp. 

35 1 24-130, June 1982 pour les transmissions numeriques et completee par 
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D.M. Duglos et R.A. Scholtz dans I'article intitule "Acquisition of spread 
spectrum signals by an adaptive array", paru dans IEEE Trans-Acous- 
Speech-Signal-Proc, vol. ASSP-37, n°8, pp. 1253-1270, August 1989. Ces 
techniques apparaissent toutefois beaucoup plus cotiteuses que les techni- 
5 ques mono-capteur et sont malheureusement inapplicables telles quelles en 
presence des derives Doppler du signal inherentes aux transmissions 
ionospheriques. 

Le but de I'invention est de pallier les inconvenients precites. 
A cet effet, I'invention a pour objet un procede permettant a un 
10 modem de se synchroniser sur un transmetteur de signaux numeriques par 
voie hertzienne en presence de brouilleurs consistant a : 

- estimer la matrice de correlation Rxx de signaux recus sur un en- 
semble de N capteurs de reception par blocs composes d'un nombre 
d'echantillons d'une sequence de synchronisation transmise par le transmet- 

15 teur, 

- a calculer la matrice de correlation inverse R^, 

- a calculer des vecteurs d'intercorrelation r Xd entre les signaux X(k) 
recus sur I'ensemble des N capteurs et un signal de replique d(k) connu, 

- a calculer un critere de synchronisation effectuant le produit scalaire 
20 entre le vecteur d'intercorrelation r Xd et le vecteur obtenu par produit de 

I'inverse de la matrice de correlation R^ et du vecteur d'intercorrelation r Xd 

- a comparer la valeur du critere obtenu a une valeur de seuil rj 
determinee pour placer la synchronisation sur I'echantillon pour lequel la 
valeur des criteres depasse la valeur de seuil r\, caracterise en ce qu'il con- 

25 siste pour estimer le decalage Doppler inherent aux communications entre 
emetteurs et recepteurs, a decouper la sequence de synchronisation en un 
nombre determine de M tranches de k 0 symboles, a effectuer sur ces tran- 
ches des calculs de vecteurs d'intercorrelation partielles r Xdi (i = 1,...,M) de 
telle sorte que sur chacune des tranches de k 0 symboles, le decalage 

30 Doppler fasse tourner la phase d'une valeur inferieure a 90°, a effectuer des 
calculs de transformee de Fourier rapide (FFT) sur les vecteurs f Xdi obte- 
nus, et a calculer une valeur du critere de synchronisation a partir des vec- 
teurs r Xd (kAf), obtenus en resultat des calculs de transformee de Fourier 
rapide, et de la matrice de correlation R^. 
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L'invention a egalement pour objet un dispositif pour la mise en 
oeuvre du procede. 

Le procede et le dispositif selon l'invention ont pour avantages qu'ils 
permettent a la fois d'assurer la synchronisation d'un modem dans des 
5 conditions de brouillage severe et en presence de multi-trajets, d'estimer la 
derive en frequence entre emetteur et recepteur et de detecter les instants 
d'arrivee de tous les multi-trajets associes au signal utile, pourvu que ceux-ci 
soient de puissance suffisamment forte par rapport a la puissance de bruit 
de fond. 

10 D'autres caracteristiques et avantages de l'invention apparaitront 

dans la description qui suit faite au regard des figures des dessins annexes 
qui represented : 

- La figure 1 un schema fonctionnel d'une antenne adaptative selon 
l'invention. 

15 - La figure 2 un exemple de mise en oeuvre d'un filtre adapte spatial a 

discrimination par la modulation pour la realisation d'une antenne adaptative 
selon l'invention. 

- La figure 3 un schema d'antenne adaptative pour la mise en oeuvre 
du procede de synchronisation spatiale selon l'invention. 

20 - La figure 4 un graphe pour illustrer la methode d'estimation de la 

matrice de correlation de signaux provenant de N capteurs. 

- La figure 5 un graphe pour illustrer le decoupage des symboles de 
la sequence de synchronisation par tranche pour le calcul de vecteur d'in- 
tercorrelation en presence de decalage Doppler. 

25 - La figure 6 un graphe montrant revolution de critere de synchroni- 

sation mis en oeuvre par l'invention en presence d'un decalage Doppler. 

- La figure 7 un organigramme retracant I'ensemble des etapes ne- 
cessaires a la mise en oeuvre du procede de synchronisation selon l'inven- 
tion. 

30 Le procede de synchronisation mis en oeuvre par l'invention fait 

appel aux techniques connues de synchronisation mono-capteur pour les 
appliquer a la realisation d'antennes adaptatives sur replique. 

Selon la technique de synchronisation mono-capteur, le critere de 
synchronisation C(n) est obtenu en calculant la correlation y(n) entre K sym- 

35 boles connus dans une sequence de synchronisation et les K echantillons 
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qui suivent I'arrivee d'un echantillon n sur le capteur, en normalisant les 
resultats par la puissance p(n) du signal recu. Ces calculs sont effectues en 
appliquant les relations : 

y(n) = ±£x(k+n)d*(k) (1) 
K k=i 
1 K 

5 p(n) = -iXx(k + n)x*(k + n) (2) 

K k=l 




ou : - x(k) est le signal issu du capteur utilise, 

- d(k) est le signal de replique connu (K symboles places en debut de 
trame), 

10 - * represente I'operation de conjugaison sur les nombres complexes. 

Ces calculs sont reiteres pour chaque echantillon n, et lorsque le cri- 
tere C(n) depasse un certain seuil r\, dependant de la probability de fausse 
alarme, la synchronisation est decidee. L'instant de synchronisation est 
determine par le rang d'echantillon n 0 pour lequel C(no)>r|. La normalisa- 

15 tion par la puissance p(n) permet de choisir un seuil independant du niveau 
de signal recu. 

Dans ce procede la synchronisation est obtenue lorsque le rapport 
signal a bruit est suffisant. Ses performances dependent en partie du nom- 
bre d'echantillons pris en compte pour la correlation. Cependant en pre- 

20 sence de brouilleurs forts, la synchronisation ne peut etre effectuee. 

L'idee mise en oeuvre dans I'invention consiste alors a effectuer la 
synchronisation, en sortie d'un traitement d'antenne adaptative, permettant 
la rejection des brouilleurs. 

Le schema fonctionnel d'une antenne adaptative est represente a la 

25 figure 1, celle-ci comporte un reseau 1 de N capteurs ^■^ ... 1 N , un reseau de 
filtres adaptes 2 (constitues par des circuits de ponderations amplitude- 
phase, des filtres a reponse impulsionnelle finie RIF ou infinie Rll) et par un 
processeur de commande adaptatif 3. Le processeur 3 calcule les coeffi- 
cients a W N du reseau de filtres 2 pour I'adapter au signal utile et aux 

30 brouilleurs, en I'absence de toute connaissance a priori sur les brouilleurs et 
a partir d'une information minimale sur le signal utile. Les coefficients a 
W N permettent de former un trou dans le diagramme de rayonnement du 
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reseau d'antenne dans la direction des brouilleurs, tout en maintenant un 
certain gain dans la direction du signal utile. 

Differents types d'information sur le signal utile sont exploites par 
I'antenne adaptative. Ceux-ci sont d'ordre spatial lorsque la direction d'arri- 
vee du signal utile est connue, d'ordre temporelle lorsque les instants de 
presence et d'absence de signal utile sont connus par exemple en evasion 
de frequence, relatifs a la forme d'onde lorsqu'une information sur la modu- 
lation du signal utile est disponible pour disposer d'une replique du signal 
utile en intercalant par exemple des sequences connues dans la trame utile, 
d'ordre spectral lorsque le signal utile occupe une certaine bande de fre- 
quence et que les brouilleurs occupent une bande plus importante, ou sta- 
tistiques pour separer les differentes sources arrivant sur le reseau en 
exploitant par exemple les statistiques d'ordre 4 des signaux en supposant 
que les differents signaux sont statistiquement independants. 

Dans le cas de I'antenne adaptative representee a la figure 2, I' infor- 
mation sur le signal utile est constitute par les sequences de symboles 
connus intercalees dans la trame entre les sequences de symboles d'infor- 
mation. Les filtres places en sortie de chacun des capteurs 1 1 a 1 N sont 
formes respectivement par des circuits de ponderation amplitude-phase 
a W N des signaux x^t) a x n (t) fournis par les capteurs. Les signaux obtenus 
en sortie des circuits de ponderation sont appliques a un circuit de somma- 
tion 4. 

Les caracteristiques du reseau de filtres 2 sont calculees en minimi- 
sant une Erreur Quadratique Moyenne (EQM) estimee £ entre un signal de 
replique d(k) et le signal y(k) sortant de I'antenne de la fagon schematisee a 
la figure 2 ou les elements homologues a ceux de la figure 1 sont reperes 
avec les memes references, i est definie par la relation 



ou : - y(k) est la sortie de I'antenne adaptative et est definie par 



- W est le vecteur poids definissant I'antenne adaptative, W* etant le 
vecteur transpose conjugue correspondant. (le signe t en exposant definit 
I'operateur correspondant.) 

- X(k) est le vecteur forme par les signaux re?us sur les capteurs. 




(4) 



y(k) = W t X(k) (5) 
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- d(k) est le signal de replique connu. 

La difference y(k) - d(k) est obtenue sur la figure 2 a la sortie d'un cir- 
cuit soustracteur 5. 

Le vecteur poids W est calcule par la formule connue de Wiener : 
5 W = R^r Xd (6) 

avec : 

1 K 

Rxx = 77 £ x M- xt M (7) matrice de correlation des signaux pro- 
K k=i 

venant des capteurs 1 1 a 1 H . (X*(k) etant le vecteur transpose conjugue du 
vecteur X(k)). 

1 K 

10 0(d = _ Z X M d M vecteur d'intercorrelation entre X(t) et d(t) (8). 

K k=i 

L'antenne adaptative ainsi utilisee conduit a une maximisation du 
rapport signal/(bruit + interference) connu sous la designation SNIR abre- 
viation anglo-saxonne de "Signal to Noise + Interference Ratio" obtenu en 
sortie de l'antenne, avec une possibility de rejection pouvant aller jusqu'a N- 

15 1 brouilleurs, et est appelee Filtre Adapte Spatial (FAS). 

La recherche de la synchronisation est effectuee selon I' invention en 
sortie de l'antenne adaptative en discriminant le signal selon sa modulation, 
et en utilisant le critere de synchronisation mono-capteur. Le schema de 
principe de la synchronisation spatiale est represente a la figure 3 ou les 

20 elements homologues a ceux des figures 1 et 2 sont represents avec les 
memes references. Le critere de synchronisation devient, en utilisant les 
notations precedentes : 

C(n) = ±2]y(k+n)d*(k) (9) 
K k=i 

y(k+n) represente la sortie de l'antenne adaptative conduisant a une 

25 minimisation de I'erreur quadratique moyenne £(n) estimee sur K echan- 

tillons entre la replique d(k) et le signal arhvant pris a partir de I'echantillon 

n, suivant la relation 

K 2 i 

^(n)= £|y(k + n)-d(k)| avec y(k + n) = W(n) l .X(k + n) (10) 

k=1 

avec: W(n) = Rxx" 1 (n).r Xd (n) (11) 
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Les differentes correlations sont calculees sur les K echantillons qui 
suivent I'arrivee d'un echantillon n suivant les relations : 



K 



Rxx(n) = 7lX(k + n)X t (k+n) (12) 
K k=i 

p Xd(n) = ^rSX(k + n)d*(k) (13) 
K k=l 

Lorsque I'echantillon n correspond a la position de synchronisation, 
les poids W(n) definissent le filtre adapte, et y(k+n) contient alors le signal 
utile debarrasse de la contribution des brouilleurs. Dans le cas contraire, 
aucun signal n'est correle avec la replique d(k) et y(k+n) ne contient plus 
alors que la contribution du bruit de fond. 

L'expression du critere de synchronisation spatiale se simplifie en 
ecrivant : 



C(n)4xw(n) t X(k + n)d*(k) = W(n) t 
K k=1 



±£x(k + n)d*(k)j (14) 



soit : 

C(n) = r< d (n)R xx (nr 1 r Xd (n) (15) 

De meme que pour le cas mono-capteur, ce calcul est reitere pour 
chaque echantillon n, et lorsque le critere C(n) depasse un certain seuil r\, 
dependant de la probability de fausse alarme, la synchronisation est deci- 
dee (la normalisation par la puissance est ici obtenue par I'intermediaire de 

R xx)- 

L'instant de synchronisation est alors determine par le rang de 
I'echantillon pour lequel C(n 0 ) est superieur a une valeur de seuil determi- 
nee r\. 

Les calculs correspondant sont effectues sur la figure 3 par le bloc 7 
de calcul du critere C(n) et par le detecteur de seuil 8, a partir de la 
sequence de synchronisation d(t) et du signal y(t) correles par le 
correlateur 6. 

Cela necessite pour chaque echantillon, d'une part de calculer N 
correlations entre le signal de replique d(k) et les signaux x(k) recus par les 
N capteurs (soit N fois les calculs effectues pour la synchronisation mono- 
capteur), et d'autre part, d'estimer puis d'inverser la matrice d'intercorrelation 
des signaux capteurs. 
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Aii cours de ces calculs la synchronisation et le calcul du filtre adapte 
2 sont effectues conjointement. Lorsque la valeur du critere depasse le seuil 
ti, calcule, la synchronisation est decidee et le vecteur poids definissant le 
filtre 2 s'en deduit facilement par la formule connue de Wiener. Le vecteur 
poids W peut ensuite etre utilise pour filtrer les echantillons contenant I'in- 
formation. 

En supposant que la replique d(t) est de puissance 1, les correlations 
tendent vers les valeurs suivantes, lorsque le nombre d'echantillons K 
devient tres grand : 

0(d = >/*s~ s 
Rxx = R + tc s S S l 

ou : - R est la matrice de correlation des interferences et du bruit de fond, 

- S est la matrice du vecteur directeur du signal utile, 

- k s est la puissance du signal utile recu sur un capteur. 

Le critere de synchronisation qui converge alors vers une valeur C 0 
s'ecrit encore : 

r TCgSW^S _ SNIRopt , * 
0 = 1 + *Sfe- 1 S = 1 + SNIRopt U6j 
ou SNIRopt est le SNIR asymptotique obtenu en sortie du filtre spatial 
adapte 2. 

Lorsque le SNIR est fort, le critere tend vers 1, lorsqu'il est faible, il 
tend vers 0. L'importance du choix du seuil ti, fixe la probability de detection 
de la synchronisation. Le seuil depend naturellement de la probability de 
fausse alarme souhaitee. 

Le calcul du critere de synchronisation C, de f Xd et de Rxx est effec- 
tue pour chaque position de la synchronisation et notamment la matrice Rxx 
est estimee et inversee egalement pour chaque position de la synchronisa- 
tion, la puissance de calcul demandee est elevee. Mais il est possible de 
reduire la puissance de calcul en estimant la matrice de correlation Rxx 
moins souvent sur un intervalle de temps plus important (K' echantillons au 
lieu de K avec K'>K de la fa?on representee a la figure 4) et en utilisant cette 
estimation de Rxx sur toutes les positions de synchronisation incluses dans 
cet intervalle de temps. 

Dans ces conditions la matrice Rxx s'ecrit : 
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ZBMBHtk) 
k=1 



k=K 0 



(17) 



ou : - B(t k ) est la contribution du bruit et des brouilleurs sur X(t k ), 

- S(t k ) est la contribution du signal utile sur X(t k ), 

- Kq est I'instant d'apparition du signal utile sur le palier, 

- e est un terme prenant en compte I'intercorrelation entre le signal 
B(t) et le signal utile S(t), qui tend vers 0 lorsque K devient grand. 

Pour la position de synchronisation testee : 



rxd = 



k 0 +k 

I 

.k=K 0 



X(t k )d» 



(18) 



Dans les conditions asymptotiques (avec K»1), lorsque le signal 
utile est synchronise avec la replique d(t), Rxx et r Xd se simplifient et 
s'ecrivent : 

R xx5 R + K '- K o + 1 ns SS t 



K' 



(19) 

r Xd = ^ s ( 2 °) 

Le critere obtenu en estimant la matrice de correlation Rxx sur K' 
echantillons possede comme valeur asymptotique, apres un calcul utilisant 
le lemme d'inversion matricielle : . 

SNIRopt 



C'o = 



1 + K'-K 0+ 1 SN|Ropt 



(21) 



Comme K ^ + 1 < 1, alors C 0 



> Cq 



La valeur du critere apparait dans ces conditions augmentee d'un 
facteur dependant de I'instant d'apparition du signal utile sur la frequence 
consideree, ce qui entraine une amelioration de la prise de synchronisation, 
au prix d'une puissance de calcul reduite. 

D'autre part, le nombre K' ne doit pas etre choisi trop grand, sa valeur 
est determinee pour que sur K' echantillons, I'environnement de brouillage 
puisse etre considere comme stationnaire. 

Le seuil a prendre en consideration depend de la probability de 
fausse alarme (Pfa). La (Pfa) est la probability qu'un bruit soit detecte 
comme une synchronisation. Le calcul de la valeur du seuil rj a lieu en cal- 
culant Pfa=P(C>Ti) lorsque I'antenne recoit uniquement du bruit. 
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En presence de bruit, le critere de synchronisation s'ecrit : 



H^lBk^ll^lBKBkM (22) 



t 1 



K 

En supposant que B k est un bruit gaussien complexe de puissance 
a 2 , independant entre les capteurs, lorsque K est suffisamment grand, le 
5 terme du milieu de la relation (22) devient : 



£ B k B k l = a +2 Id (23) (Id etant la matrice identite ou uni- 



taire) 

le critere de synchronisation peut alors etre approxime de la facon 
suivante : 



1 



10 c 'Ml B ^M^ 

Soit: 



(24) 



N 



K 

ZBk'dk' 
k=1 



1 



ifff + I?) (25) 



avec: £ B k \ = Rj + j L (26) 
k=1 

B k represente le bruit obtenu sur le capteur i, pour I'echantillon k. 
R, + j Ij peut etre approxime par un bruit gaussien complexe de 
variance K a 2 . Dans ces conditions Rj et Ij sont chacun des bruits gaussiens 



reels de variance . En notant T; = Rj ou h, C devient : 

2 iii 

2N, 



ou Tj est une variable aleatoire qui suit une loi normale centree de variance 



J2 



20 La variable aleatoire Zj = Tj suit done une loi normale centree 

2 VKa 2 

de variance 1 . 

C devient : 
. 2N 

C4IZ? (28) 
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Comme la suite des Z\ peut etre consideree comme une suite de 2N 
variables aleatoires independantes gaussiennes de moyennes nulles et de 

2N 

variances 1, la variable aleatoire £Z 2 suit une loi du x 2 a 2N degres de 

i=1 

liberte. 

La probability d'avoir C > n s'ecrit alors : 



p(C> n ) = p|zZ 2 >2KT 1 j = Pfa (29) 



Le tableau suivant donne differentes valeurs de seuils ti obtenues par 
recherche dans une table donnant P(Z>A) pour une loi du x 2 a 2N degres de 
liberte et pour Pfa = 10" 5 . 

N = 5 => A = 42 ti = 0,26 

N = 4 => A = 37 => ti = 0,23 

N = 3 => A = 33 ^ n = 0,21 

N =2 ^ A = 28 ^> ti = 0,18 

Comme entre remission et la reception, le signal transmis peut etre 
affecte d'un certain decalage en frequence, ou decalage Doppler, du essen- 
tiellement aux derives en frequences des oscillateurs locaux des syntheti- 
seurs a remission et a la reception, ainsi qu'aux conditions de propagation 
ionospherique, il est necessaire de compenser ce decalage. Le risque est 
s'il n'est pas compense, d'empecher toute prise de synchronisation. 

En appelant Af 0 , le decalage Doppler entre le signal recu par I'an- 
tenne et le signal attendu et en considerant la position de synchronisation, le 
vecteur signal recu peut s'ecrire avec les notations des paragraphes sui- 
vants : 

La matrice de correlation estimee peut alors s'ecrire : 
RXX =^ l[v^d(tk)e 2|7t Afo tk S + B(t k )J^d(t k )e 2 -j-' l Af o t K.s + B(t k )] t (30) 



k=r 

K 

Rxx=e+ ' 



^Z^sWtKfsS^BO^B^) 1 ] (31) 

ou s est une matrice prenant en compte les intercorrelations entre la contri- 
bution du signal utile et le bruit b(t). Lorsque K»1, dans les conditions 
asymptotiques, Rxx peut s'approximer par : 
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Rxx = n s S S + + a 2 Id (32) 

La relation (32) montre que le decalage Doppler ne modifie pas le 
calcul de Rxx- En effet, le decalage s'elimine naturellement lorsque le signal 
utile est multiplie par sa quantite conjuguee. 

Le vecteur d'intercorrelation peut s'ecrire : 

*d = r\ Z[v^ *kh 21 * ^ S + B(t k )] d (t k ) (33) 



p Xd = 6' + r\ l[v^|d(t k )| 2 e 2 ^^ s] (34) 



ou s' est un terme prenant en compte I'intercorrelation entre le bruit B(t) et la 
replique. Lorsque K»1, r X( i peut s'approximer par : 



^lV^KkT)| 2 e 



2\n. 



S (35) 



La relation (34) montre que le decalage Doppler modifie I'estimation 
du vecteur d'intercorrelation r Xd et que du fait que la phase du terme en ex- 
ponentiel complexe peut tourner de plusieurs fois 360° sur I'intervalle de 
temps utilise pour le calcul de la correlation, integration ne peut plus avoir 
lieu de facon coherente ce qui empeche le critere de synchronisation multi- 
capteurs de converger. Le processus suivant permet d'estimer le decalage 
Doppler et de le compenser pour rendre possible la prise de synchronisa- 
tion. 

Selon ce processus, les K symboles de la sequence de synchronisa- 
tion sont decoupes de la fa?on representee a la figure 5 en M tranches de 
Kq symboles (M Kq = K), de telle sorte que sur chacune des tranches de Kq 
symboles, le decalage Doppler fasse tourner la phase du signal recu par 
I'antenne d'une valeur inferieure a 90°, afin de permettre une integration 
coherente de r X( j sur chacune des tranches. 

Si r X( jj designe la valeur de I'intercorrelation estimee sur la tranche i, 
celle-ci s'ecrit : 



r Xdi ■ 



(i+i)K 0 

I 

k=iK 0 +1 



v^|d(kT)i 2 



2j*.^.k 



S i = 0,...,M-1 (36) 
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En mettant en facteur dans la relation (36) un terme dependant uni- 
quement de Kg, la longueur de la tranche elementaire consideree, et du 
decalage Doppler, la relation (36) devient : 



r X di = 



4 K 0 , 2jiiAk 

lX^|d((k + iK 0 )T)| 2 e F e 
* k=1 



e F * .S i = 0 M-1 (37) 



Ce qui peut encore s'ecrire : 
2j,AiK 0 

rxdjSXje Fe S i = 0,...,M-1 (38) 

La relation (38) montre que les vecteurs r Xdi sont des vecteurs, 

echantillonnes a Leur transformee de Fourier discrete est alors : 
K 0 

M-1 -2jn.—^i 

P X d(Af)=Ir Xdi e (39) 
i=0 

M-1 2 j*.^.iK 0 
rxd(Af)= I xje ^ S (40) 
i=0 

La relation (40) montre que lorsque Af = Af 0 , le decalage Doppler est 
compense, ce qui rend maximum la valeur du critere de synchronisation 
rxd(Af) t Rxx" 1 r Xd (Af). 

Afin d'estimer Af, c'est-a-dire afin d'estimer le decalage Doppler qui 
conduit a la valeur maximale du critere, et afin de reduire le cout de calcul 
correspondant, une transformee de Fourier rapide est effectuee sur les 
vecteurs r Xdi (le nombre de tranches elementaires M etant choisi pour former 
une puissance de 2) : ce calcul donne : 

M-1 - 2 \n.-^-.\k M-1 _2i7r — 

fa(ktf,) = I P Xd i e = I r Xdi e M k = 0 M-1 (41) 

i=0 i=0 
ou Af-] est la resolution Doppler telle que : 

Af, = ^> (42) 

Si la resolution est insuffisante, une amelioration d'un facteur 2 peut 
etre obtenue en calculant une transformee de Fourier rapide, sur les vec- 
teurs r Xdi "decalees" en frequence de la facon suivante : 
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( aM M_1 I -2m- 
r Xd kAfl^=Z p x di e M e J M k = 0 ,...,M-1 (43) 

i=0 V 

L'operation consiste a multiplier chacun des vecteurs d'intercorrela- 
tion r Xc ji par le terme e~^ l/M ^ et a effectuer ensuite un calcul de transformee 
rapide sur les vecteurs ainsi obtenus. 

Le critere de synchronisation est ensuite estime sur chacune des 
cases de calcul des transformees de Fourier precedentes par la relation : 

c ( k fH( k f) R ^ f 4f) k=o M - 1 (44) 

Bien evidemment, si la resolution est toujours insuffisante, il est pos- 
sible de I'ameliorer en utilisant le meme processus pour calculer des vec- 
teurs Pxd^kAf-i + j oil n est un entier quelconque. 

Un test sur la case de la transformee correspondant au maximum du 
critere de synchronisation est effectue. Si le seuil est depasse pour une 
valeur k 0 ^ donnee du decalage Doppler, il est alors possible d'affiner le 

calcul du decalage Doppler (indispensable pour la demodulation) en calcu- 
lant la valeur du critere pour differentes valeurs de decalage Doppler entou- 
rant k 0 . En effet, a Tissue dexette etape, la resolution sur ('estimation 
Afi 

du decalage Doppler est Pour I'ameliorer, il faut proceder par dichoto- 

mie : a chaque iteration et de la fagon representee a la figure 6, la valeur du 
critere pour la valeur de decalage Doppler situee au milieu des deux valeurs 
de decalage Doppler donnant les valeurs maxima du critere lors de I" iteration 
precedente est calculee. La recherche par dichotomie est arretee des que la 
resolution souhaitee sur I'estimation du decalage Doppler est obtenue. 
Apres n iterations de I'algorithme de recherche par dichotomie, la resolution, 

At 

qui etait initialement de — L ,est: 

i) 

L'enchainement des etapes de calcul precedente pour la mise en 
oeuvre du procede selon I'invention est figure sur I'organigramme de la 
figure 7. 
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II consiste suivant I'etape 9 a estimer la matrice de correlation Rxx 
suivant la relation (7). Cette matrice est estimee sur un bloc contenant un 
nombre d'echantillons K' superieur au nombre d'echantillons de la sequence 
de synchronisation. 

5 L'estimation de Rxx est utilisee ensuite pour toutes les positions de 

synchronisation incluses dans le bloc. Puis la matrice de correlation inverse 
est calculee a I'etape 10. r Xdi qui est defini par la relation (8) est calcule 

a I'etape 1 1 . Pour chaque echantillon correspondent a une position de syn- 
chronisation P dans la suite des echantillons P = 1 a P = K'-K, r Xd j est calcu- 
10 le pour i = 1 a M sur M blocs de Kq symboles comme defini a la figure 5 sui- 
vant la relation 

p Xdi = it- S X(P + K + iKo) d (K + iKo) (46) 
K 0 K=1 

Des calculs de transformee de Fourier Rapide sont effectues a I'etape 
suivante 12 sur un nombre determine M (M = 8 par exemple) de vecteurs 
15 r X( jj suivant la relation (41 ). 

Aux etapes 13 et 14 un calcul de vecteur r X{ jj decale et un calcul de 
transformee de Fourier rapide suivant la relation (43) sont executes pour 
ameliorer la resolution par un facteur de 2. 

Le critere de synchronisation est calcule a I'etape 15 en application 
20 de la relation (44). 

La comparaison du critere ou seuil r\ est effectuee a I'etape 16. Si a 
I'etape 16 la valeur du critere c ^~^j est superieure ou egale au seuil, la 

position de synchronisation est definie a I'etape 17 sur la position P de 
I'echantillon correspondant. L'estimation du decalage Doppler peut ensuite 
25 etre affine a I'etape 18 en calculant a nouveau la valeur du critere de syn- 
chronisation pour des decalages Doppler entourant le decalage ^k 0 

precedent. Cette etape peut etre suivie par une etape 19 d'estimation de la 
position de synchronisation pour les autres trajets. 

La recherche de synchronisation precedemment decrite est effectuee 
30 pour chacun des echantillons de signal numerise. En pratique, pour amelio- 
rer la precision de la position de la synchronisation il est preferable d'effec- 
tuer un surechantillonnage du signal d'un facteur determine L, L etant egal a 
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4 par exemple. Les traitements de I'organigramme de la figure 7 et decrits ci- 
dessus sont encore utilisables sur le signal ainsi surechantillonne, mais il est 
possible d'utiliser une variante dans le procede, afin de reduire la puissance 
de calcul necessaire au traitement. Dans ces conditions, apres avoir calcule 
5 pour chaque echantillon le vecteur r Xdj pour i = 1, ... M sur M blocs de Kq 
symboles suivant la relation (46) un critere de synchronisation Cj est calcule 
sur chacun des M blocs par la relation 
Ci = r< d ,R^.r Xdi (47) 

Puis un critere somme de I'ensemble des Cj est calcule suivant la 
10 relation : 

Ci(P) = £Cj (48) 
i=1 

Ce critere correspond a une integration non coherente des Ci 
(contrairement a ce qui est fait pour le calcul par FFT ou on integre de facon 
coherente), ce qui se traduit par une remontee du critere sur les echantillons 
15 de bruit. 

Le calcul du critere a lieu par blocs de L echantillons consecutifs. Sur 
chacun de ces blocs, I'echantillon qui donne la valeur maximum de Ci(P) est 
selectionne pour effectuer le calcul par transformee de Fourier rapide 
necessaire a I'estimation du decalage Doppler comme deja decrit aux etapes 
20 1 1 a 1 9 de la figure 7. 
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REVENDICATIONS 

1 . Procede permettant a un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de 

5 brouilleurs consistant a : 

- estimer (9) la matrice de correlation R^ de signaux recus sur un 
ensemble de N capteurs de reception (1-j ... 1 N ) par blocs composes d'un 
nombre d'echantillons d'une sequence de synchronisation transmise par le 
transmetteur, 

10 - a calculer (10) la matrice de correlation inverse R^, 

- a calculer (11) des vecteurs d'intercorrelation f Xd entre les signaux 
X(k) recus sur I'ensemble des N capteurs et un signal de replique d(k) 
connu, 

- a calculer un critere de synchronisation (7) effectuant le produit 
15 scalaire entre le vecteur d'intercorrelation r Xd et le vecteur obtenu par pro- 
duit de I'inverse de la matrice de correlation et du vecteur d'intercorre- 
lation r Xd 

- a comparer (8) la valeur du critere obtenu a une valeur de seuil r] 
determinee pour placer la synchronisation sur I'echantillon pour lequel la 

20 valeur des criteres depasse la valeur de seuil ri, caracterise en ce qu'il con- 
siste pour estimer le decalage Doppler inherent aux communications entre 
emetteurs et recepteurs, a decouper la sequence de synchronisation en un 
nombre determine de M tranches de k 0 symboles, a effectuer (11, 13) sur 
ces tranches des calculs de vecteurs d'intercorrelation partielles 

25 rxdi (i = 1, .... M) de telle sorte que sur chacune des tranches de k 0 symboles, 

le decalage Doppler fasse tourner la phase d'une valeur inferieure a 90°, a 
effectuer (12, 14) des calculs de transformee de Fourier rapide (FFT) sur les 
vecteurs r Xdi obtenus, et a calculer une valeur du critere de synchronisation 
a partir des vecteurs r Xd (kAf), obtenus en resultat des calculs de 
30 transformee de Fourier rapide, et de la matrice de correlation R^ . 

2. Procede selon la revendication 1 caracterise en ce que I'estimation 
de la matrice de correlation R^x est obtenue par la relation 

1 K 

R)0( =-^x(k).x^) et le vecteur d'intercorrelation est estime suivant la rela- 
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1 K 

tion r X( j = — X^k+p)^)' ou p est la P os ' tion du Premier echantillon a partir 

duquel la correlation est effectuee entre le vecteur X et la replique d. 

3. Procede selon la revendication 2 caracterise en ce qu'il consiste a 
choisir la case de calcul de la transformee rapide qui donne une valeur 

5 maximum du critere de synchronisation et a comparer la valeur maximum 
obtenue a une valeur de seuil determinee pour decider de la synchronisa- 
tion. 

4. Procede selon Tune quelconque des revendications 1 a 3 caracte- 
rise en ce qu'il consiste pour obtenir une meilleure resolution sur le deca- 

10 lage Doppler et une valeur plus importante du maximum du critere, a calcu- 
ler (13) des vecteurs r x< jj decales en frequence et a effectuer un calcul de 
FFT sur les vecteurs r Xdi decales de facon a obtenir la valeur du critere avec 
une meilleure resolution sur les frequences Doppler. 

5. Procede selon Tune quelconque des revendications 1 a 4 caracte- 
15 rise en ce qu'il consiste a surechantillonner les signaux recus sur les N 

capteurs par rapport au rythme symbole, pour ameliorer la resolution sur la 
position de synchronisation. 

6. Procede selon Tune quelconque des revendications 1 a 5 caracte- 
rise en ce qu'il consiste a effectuer le calcul de la matrice de correlation Rxx 

20 sur un nombre d'echantillons K' superieur au nombre K d'echantillons de la 
sequence de synchronisation afin de diminuer la puissance de calcul. 

7. Procede selon I'une quelconque des revendications 1 a 6 caracte- 
rise en ce qu'il consiste pour diminuer la puissance de calcul et apres avoir 
surechantillonne d'un facteur L par rapport au rythme symbole, a effectuer 

25 une etape de pretraitement consistant a calculer pour chaque echantillon p 
(p = 1 a L) composant un symbole les valeurs Cj du critere de synchronisa- 
tion sur chacune des M tranches de k 0 symboles et a calculer la valeur du 
critere en effectuant la somme des criteres Cj obtenus pour retenir dans la 
suite du traitement I'echantillon donnant la valeur maximale du critere 

30 C = XCj. 
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Estimation de R xx sur K' echantillons 
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